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In this work, we have considered emergent universe in generalized gravity theories like Chameleon, 
f{R) and f{T) gravities. We have reconstructed the potential of Chameleon field under emergent 
scenario of the universe and observed its increasing nature with evolution of the universe. We 
have revealed that in the emergent universe scenario, the equation of state parameter behaves like 
quintessence in the case of f{R) gravity and like phantom in the case of / (T) gravity. 

PACS numbers: 

I. INTRODUCTION 

The central paradigm for structure formation in the universe is the inflationary scenario. Under very general 
conditions, inflation is future eternal, in the sense that once inflation has started, most of the volume of the 
universe will remain in an inflating state [1]. We consider the total density U,{t) = 1 + ^^j^^ , where, k, a and 
H denote the curvature, scale factor and Hubble parameter respectively. Inflation drives the curvature term 
towards 0, which docs not imply that fc = 0. However, the standard inflationary model is based on a flat 
{k = <^=^ ^0 — 1) Friedmann-Robertson- Walker (FRW) geometry, motivated by the fact that inflationary ex- 
pansion rapidly wipes out any original spatial curvature. Leading alternatives to inflation that are motivated by 
recent advances in string/M-theory are the pre-big bang [2] and ekpyrotic/cyclic [3] scenarios, respectively. The 
spatial curvature of the real universe is in principle determined by observations. The recent WMAP data seems 
to point to a universe that is close to (but not quite) flat, with a total density parameter of fttot = 1-02 ± 0.02. 
Theory will have to give way to data if the data clearly tell us that ilo > 1 ([4], [5]). If JIq is taken as 1.02, then 
the power spectra of CMB anisotropics and matter can show testable differences from the standard flat model [6] . 

The singularity theorems assume that either the universe has open space sections, implying /c = or — 1; 
or the Hubble expansion rate H = ^ is bounded away from zero in the past. There are inflationary universes 
that evade these constraints and hence avoid the conclusions of the theorems [7]. Ellis and Maartens [4] 
ITJi considered closed models in which k = +1 and H can become zero, so that both of the assumptions of the 
inflationary singularity theorems are violated. In this paper [4], it was shown that if A; = +1 then there are 
I closed inflationary models that do not bounce, but inflate from a static beginning, and then reheat in the usual 
way. 



The inflationary universe emerges from a small static state that has within it the seeds for the development 
of the microscopic universe and it is called Emergent Universe scenario. The universe has a finite initial size, 
with a finite amount of inflation occurring over an inflnite time in the past and with inflation then coming to 
an end via reheating in the standard way. Reference [8] (see also [11]) summarized the features of emergent 
universe as: 

1. the universe is almost static at the finite past (t —oo) and isotropic, homogeneous at large scales; 

2. it is ever existing and there is no timelike singularity; 

3. the universe is always large enough so that the classical description of space-time is adequate; 

4. the universe may contain exotic matter so that the energy conditions may be violated; 
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5. the universe is accelerating as suggested by recent measurements of distances of high redshift type la 
supernovae. 

Ellis ct al [9] provided a realization of a singularity-frcc inflationary universe in the form of a simple 
cosmological model dominated at early times by a single minimally coupled scalar field with a physically 
based potential. Mukherjee et al [8] presented a general framework for an emergent universe scenario 
and showed that emergent universe scenarios are not isolated solutions and they may occur for different 
combinations of radiation and matter. Campo ct al [10] studied the emergent universe model in the context of 
a self-interacting Jordan-Brans-Dicke theory and showed that the model presents a stable past eternal static 
solution which eventually enters a phase where the stability of this solution is broken leading to an inflationary 
period. Debnath [11] discussed the behaviour of different stages of the evolution of the emergent universe 
considering that the universe is filled with normal matter and a phantom field. Mukherji and Chakraborty [12] 
developed Einstein-Gauss-Bonnet (EGB) theory in the emergent universe scenario. They have considered the 
Friedman-Lematre-Robertson- Walker cosmological model in Horava gravity and the emergent scenario for all 
vahies of the spatial curvature. Paul et al [13] predicted the range of the permissible values for the parameters 
associated with the constraints on exotic matter needed for an emergent universe. The emergent universe in 
Horava gravity was studied by Mukherji and Chakraborty [14]. 

Many theoretical models that anyhow describing the accelerated expansion of the universe and which 
appears to fit all currently available observations are affected by significant fine-tuning problems related to 
the vacuum energy scale and therefore it is important to investigate alternatives to this description of the 
Universe [15]. There exist several other approaches to the theoretical description of the accelerated expansion 
of the universe. One of these is a modified gravity theories. Studies of the physics of these theories is however 
hampered by the complexity of the field equations, making it difficult to obtain both exact and numerical 
solutions which can be compared with observations. These problems can be reduced somewhat by using the 
theory of dynamical systems, which provides a relatively simple method for obtaining exact sohitions and a 
description of the global dynamics. Modified gravity constitutes an interesting dynamical alternative to ACDM 
cosmology in that it is also able to describe with success the current acceleration in the expansion of our 
universe, the present dark energy epoch [16]. Presumably the simplest modified gravity model of dark energy 
is so-called ,f{R) gravity in which / is a function in terms of a Ricci scalar R ([17], [18], [19], [20], [21], [22]). 
The ,f{R) theories do not seem to introduce any new type of matter and can lead to late time acceleration [23]. 
Another modified gravity theory is the Chameleon field theory where the scalar field interacts with some kind 
of matter, behaving as perfect fluid. The Lagrangian of the model contains a term, where the matter interacts 
with the effective metric (physical metric, multiplied by a conformal factor depending on the scalar field) ]24]. 
Brax et al [23] reviewed both f{R) gravity and Chameleon field theories and discussed association between 
them. Another interesting sort of modified theories is so-called f{T) gravity (where T is the torsion) ([25], 
[26], [27]). Recently, it is shown that such f(T) gravity theories also admit the accelerated expansion of the 
universe without resorting to dark energy ([27], [28]). In the present paper, we have discussed the emer- 
gent universe scenario in the generalized gravity theories like Chameleon gravity, f{R) gravity and /(T) gravity. 

The organization of the paper is as follows: In section II, we have given the basic field equations in the 
Einstein's gravity. We have chosen some particular form of scale factor for emergent scenario of the universe 
and made some restrictions upon the parameters. In section III, we have discussed the Chameleon field gravity 
model and shown the behaviour of the potential function for emergent scenario. In section IV and V, we have 
studied the ,f(R) and /(T) gravity models respectively and shown the nature of equation of states in emergent 
scenarios of the universe. Finally, some conclusions have been drawn in section VI. 



II. BASIC EQUATIONS 

For a FRW spacetime, the line element is given by 



dr"^ 



1 — fcr^ 



r^{de^ 



sm 



(1) 



where a{t) is the scale factor and k (= 0,±1) is the curvature scalar. Now consider the Hubble parameter 
{H) and the deceleration parameter (g) in terms of scale factor as 
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So the Einstein field equations are given by (choosing c = 1) 
and 



H-^=-AttG{p + p) (3) 



and the energy conservation equation is 



p + 2,H{p + p)^Q (4) 
where p and j» are respectively the density and pressure of the fluid. 

For emergent universe, the scale factor is chosen as [8] 

a = ao(/3 + e"*)" (5) 
where the constant parameters are restricted as follows [14]: 

1. ao > for the scale factor a to be positive, 

2. /3 > 0, to avoid any singularity at finite time (big-rip), 

3. a > or n > for expanding model of the universe, 

4. a < and n < implies big bang singularity sX t = —oo. 

For this choice of the scale factor, the Hubble parameter, its derivatives are given by 



~ e«* + /3 ' - (e«* + /3)2 ' ~ (e«* + 

Here, H, H, and H tend to as t ^ — oo. We have seen that i? > for emergent scenario. Prom (2) and 

(3), it can be observed that the emergent scenario is possible for (i) flat (k = 0) universe in phantom stage only, 
(ii) open (A; = — 1) universe only for phantom stage ii H > ^ and (iii) closed {k = +1) universe in quintessence 
stage i{ H < ^ and in phantom stage if iJ > ^ . 

Also we calculated the cosmographic parameters [29, 30] 



^ a at-' 



J = 



— 1 d^g 



a dt^ 



„ _ 1 d*a fJ-4 



a at^ 



(6) 



which are usually referred to as the deceleration, jerk, snap, and Icrk parameters. After calculating the 
cosmographic parameters for the emergent universe we have plotted them against cosmic time t in the figures 
la and lb. It is observed in figure la that the deceleration parameter is staying at negative level. This indicates 
the ever accelerating nature of the emergent universe. In figure lb we observe that the snap, jerk and lark 
parameters are increasing function of the cosmic time in the emergent universe. 



Fig. la 



Fig. lb 



Fig. la shows the deceleration parameter q for emergent universe. Fig. lb shows snap parameter s (green line, 
middle), jerk parameter j (red line, bottom) and lerk parameter I (blue line, top) for emergent universe. 



III. EMERGENT SCENARIO IN CHAMELEON GRAVITY 

In the flat homogeneous universe, we consider that the relevant action given by [31, 32, 33] 



S = 



-gd X 



1 



n 



(7) 



where is the Chameleon scalar field and V{(j)) is the Chameleon potential. Also, TZ is the Ricci scalar and G 
is the Newtonian constant of gravity, f{(j))C is the modified matter Lagrangian and /(</>) is an analytic function 
of <j>. The variation of action (9) with respect to the metric tensor components in a FRW cosmology yields the 
field equations (choosing 8nG =1), 



and 



H 



Pf + + v{4>) 



(8) 



(9) 



The conservation equation and the wave equation in presence of Chameleon field are given by [31] 



-{pf) + 3Hip + p)f = {p + p)f (10) 

and 

3H^'^ + <j)4> + V+{p + p)f^0 (11) 

The early universe contains a standard scalar field cj) with energy density p^ = ^(j>^ + V{(f>) and pressure = 

^(/)^ — V{(t>) and possibly also ordinary matter with energy density p and pressure p = wp, where — | < w < 1. 
The cosmological constant is absorbed into the potential V. There are no interactions between matter and the 
scalar field. Following references [31,32] we make the simplifying assumption 



V = Vo<P^ (12) 

where, Vq is a positive constant. Using this simplifying assumption in the Chameleon field equation along 
with the specified form of scale factor (5) and equations (8) and (9), we get (after simplification) the form of (j> 
in the emergent scenario as 
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Fig. 2a. This figure shows the variation of potential V of the Chameleon field with cosmic time t in the 
emergent universe scenario. The cases of flat (fc = 0) (the blue line, middle), open (fc = —1) (the red line, bottom) and 

closed {k = 1) (the green line, top) universe are considered. 
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Fig.2b Fig.2c 



Fig. 2b shows the variation of the analytic function / of the Chameleon field with cosmic time t in the emergent 
universe scenario. In this figure, the cases of flat {k — 0) (the blue line, top), open (fc = — 1) (the red line, middle) and 
closed {k = 1) (the green line, bottom) universe are considered. Fig.2c shows the variation of the analytic function / of 
the Chameleon fleld with <j) in the emergent universe scenario. In this figure, the cases of flat (fc = 0) (the blue line, 
bottom), open {k — —1) (the red line, top) and closed (fc — 1) (the green line, middle) universes are considered. 
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Fig.3 

Fig. 3 represents the evolution of the scalar field (j) of the Chameleon field with time t in the emergent universe scenario 
for open (fc = — 1) (the red line, top), closed (k = 1) (the green line, bottom) and flat (fc = 0) (the blue line, middle), 

universe respectively. 
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Fig. 4a Fig. 4b 

Fig. 4a shows the potential V against the scalar field </!> in the cases of {k — —1) (the red line, bottom) and closed 
(fc = 1) (the green line, top). Fig. 4b shows the potential V against the scalar field (f} in the case of flat k — universe. 



' 2(fc(l + 3fcw)(e°* + /?)-2" + e"*7ia2(e"* + /3)-2(3e"*n(l + w) + 2/3)ag) 

al{-l + w + 2{l + w)Vo) ^ ' 

Using (5), (8), (12) and (13) and the p = /9oa-3(i+™)/"' (from (10)) for Chameleon field from reference [31], 
we express /(</>) as a function of t as follows 



/ 



((e*" + /3)"ao)i+3"'(fc(e*" + /3)2(-4 + 3(1 - k)w + (4 + 6(1 - k)w)Vo) - 2e*°'na^{e*°' + /3)2"ag(3e*"n + (3(1 + 2Va))y 



((e*" + /3)2(-l + w + 2{l + w)Vo)po) 

(14) 

We see that (e"* + /3) 2(3e"*n(l + w) + 2/3)a§) is always positive. Thus from (13) we see that emergent 
scenario is possible for flat, open and closed universe under the following conditions: 

• In the case of fc = 0, 0^ > q for Vo > 



2(l+«;) 



. In the case of fc = -1, > q for > and Vn > 



3(l+e2"tn2Q,2a2(g<,t_|,^)2(-l + „)) aii^ "U ^ 2{l+w) 

• In the case of fc = 1, 0^ > for ■(« > ~ ^~^''^^^"J^t„^j^^^^(^Jtlf)2i^i^^^^ and Vq > j^j^- 

Thus, under the following conditions we consider the cases of flat (fc = 0), open (fc = —1) and closed (fc — +1) 
universe. In figure 2a, we plot the potential against cosmic time. We find that the potential is increasing with 
time in all of the three cases. While plotting, the constant parameters are chosen according to the above four 
constraints. Evolution of f{(p) is presented in figures 2b and 2c. In all of the three situations f{(p) is staying 
at negative level and is showing a falling behaviour with cosmic time as well as the scalar field (j). For open 
(fc = —1), closed (fc = +1), and fiat (fc — 0) universe, the evolution of the scalar field </> exhibits an increasing 
pattern with evolution of the universe. These are presented in figure 3 where the scalar field (j) is plotted 
against cosmic time t. Figure 4a shows that potential V is decreasing with the scalar field cj) for fc = — 1 and 
fc = +1. However, figure 4b shows that for fc = i.e. for open universe, V is increasing with (j). 



IV. EMERGENT SCENARIO IN f{R) GRAVITY 



Recently, motivated by astrophysical data which indicate that the expansion of the universe is accelerating, 
the modified theory of gravity (or f{R) gravity) which can explain the present acceleration without introducing 
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dark energy, has received intense attention [34-40]. In a recent work, Nojiri and Odinstov [38] suggested 

two realistic f{R) and one F{G) modified gravities characterized by the presence of the effective cosmological 
constant epochs in such a way that early-time inflation and latc-time cosmic acceleration arc naturally unified 
within single model. In another work, Nojiri and Odintsov [39] proposed another class of modified f{R) gravity 
which unifies i?™ infiation with ACDM era. Cognola et al [40] demonstrated a class of viable modified f{R) 
gravities describing infiation and the onset of accelerated expansion. In f{R) gravity action is described by 
an arbitrary function of the scalar curvature R. Extensive review of f{R) gravity is available in Nojiri and 
Odintsov [21, 22]. The action of f{R) gravity is given by [44-47] 



m 

2^2 



(15) 



where g is the determinant of the metric tensor g^^, Cmatter is the matter Lagrangian and = SnG. The/(i?) 
is a non-linear function of the Ricci curvature R that incorporates corrections to the Einstcin-Hilbcrt action 
which is instead described by a linear function f{R)- The gravitational field equations in this theory are 



3/'(i?) 



(P + Pc) 



(16) 



H 



2/'(i?) 



iP + P + Pc+Pc 



(17) 



where pc and pc can be regarded as the energy density and pressure generated due to the difference of f{R) 
gravity from general relativity given by [45] 



-(-f{R) + Rf'{R))-3HRf"{R) 



(18) 



Pc= ^ 



-,ifiR) 



Rf'{R)) + {2HR + R)f"{R) + R?f"'{R) 



(19) 



where, the scalar tensor R = 6{H + 2H'^) [39,45]. Here, if = | is the Hubble parameter. Here, p and p are 
the matterenergy density and pressure respectively. In the spatially flat universe, i.e. k = 0, with the matter 
as dust, namely p = 0, following Feng [45], p = 3HQil,rnoe~^^ with x = Ina, flmO = Pmo/3-ffo is the 

present Hubble parameter. Following references [36] we choose f{R) as 



f{R) = R + £^Rt' + QR-^ 
where, constants > 0; v > Subsequently, we get from equation (17) 



(20) 



H = {a^HR{(i^{l + 1^)+ Rf'+'^i-l + fi)ijOR H 
a^R{-Cu{l + 1^)- i?''+'^(- 



-l + p)pi)R 



■SR^+'^K'^H^flr. 



-2 4 
>o)x 



-p){-i + p)pOr^ 



(21) 



{2a-^R^{-Cu + R^iR + iJ^MO))" 



As both R and H are functions of the scale; factor a and its derivatives, it seems that the equation (21) is 
set for yielding the solution for the scale factor. However, it involves fourth-order derivatives of a {R already 
contains a) and is highly nonlinear. This makes it difficult to obtain a completely analytic solution for a. 
Mukhcrjce ct al [8] obtained the general solution of the scale factor for emergent universe in the form presented 
in equation (6) without referring to the actual source of the energy density. Se we use equation (6) as the 
choice of scale factor for emergent nature of the universe and reconstruct pc and Pc based on R in the emergent 
universe and are expressed as follows 
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Pc = (2„e'.+ff).^. (2-^-"3-'-g-'-((l + i.)(-4e2t«n2 - e*"(-z. + 4n(l + v))P - (1 + z.)/32)C 

_^^M+i')(22+M+>'3A'+i'e2t"n2(-l + - e^+^e*" (22) 
X (4n(-l + /i) - + /i)/? + (2i+^+''3'^+''// - 6^+''(l + 



»c(2ne'°+ff)^ (6"^e-^"/3(e^"(-l + 4n) + /3)(e*"(l + 2na) + /3)x 
+ (1.(1 + i.)C + 6''+'^(-l + + „(2ne'^+ff)a 6"^e'^"/3'(e*°(-l + 4n) + /3)2(-i.(l + :.)(2 + ^)C+ 

+6^+''(-2 + ^)(-l + 

(23) 

where 



na2e*"(2ne*« + /3) 
(e*" + /3)2 



(24) 



The Pc and pc are plotted in figures 5a and 5b respectively. The figures show that although pc is increasing 
with cosmic time t, the negative Pc is falling with cosmic time. While summarizing the attractiveness of the 
modified gravity approach, Nojiri and Odintsov [35] stated that "modified gravity may naturally describe 
the transition from non-phantom phase to phantom one without necessity to introduce the exotic matter" . 
Abdalla et al [49] had shown in an earlier study that the modified gravity minimally coupled with the usual 
(or quintessence) matter may reproduce the quintessence (or phantom) evolution phase for the dark energy 
universe in an easier way than without such coupling. In a more recent study, Bamba et al [50] studied a 
viable model of modified gravity in which the transition from the de Sitter universe to the phantom phase 
can occur. Motivated by these earlier studies we examine the bahaviour of the equation of state parameter in 
/(i?) gravity. The equation of state parameter plotted in figure 6 is found to stay above —1. This indicates 
quintessence behaviour. This means that evolution from quintessence from phantom is not available in the case 
of f{R) gravity in emergent universe. From figure 7 we find the increasing behavior of f{R) with the evolution 
of the universe. 



V. EMERGENT SCENARIO IN /(T) GRAVITY 



Recently, models based on modified teleparallel gravity were presented as an alternative to inflationary models. 
The theory so obtained is called as the /(T) theory [51-57]. In this theory instead the curvature defined via 



Fig.( 



Fig.7 



Fig. 6 shows the evolution of the equation of state lUc = — with cosmic time t. 
Fig.7 represents the evolution of f{R) with cosmic time t. 



the Levi-Civita connection, the so-called Weitzenbock connection is used. But in this case the theory has no 
curvature but instead torsion. Similar to general relativity where the action is the curvature scalar, the action 
of teleparallel gravity is a torsion scalar T ([51], [52], [58] ). The action of /(T) theory is obtained by replacing 
T in the action of teleparallel gravity by T + f{T) [53]. We start with the following action for the /(T) gravity 
[52] 



(25) 



where T is the torsion scalar, /(T) is general differentiable function of the torsion and Cm corresponds to the 
matter Lagrangian and ~ SttG. The torsion scalar T is defined as 



— '-'p ^ pu 



with 



(26) 



(27) 



where, K'^'^ is the contorsion tensor, which equals the difference between Weitzenbock and Levi-Civita 
connections. So /(T) gravity uses the curvatureless Weitzenbock connection. 

The modified Friedmann equations of motion are ([27], [56]) 



SttG 



f{T) 



2/'(r)ij2 



(28) 



i? = - 



ATTGjpm + P,n) 

1 + f'{T) - f"{T) 



(29) 



where, Pm and Pm respectively for the energy density and pressure of the matter content of the universe, 
with equation-of-state parameter Wm = Pm/ Pm- Since we are interested in the late time universe, we consider 
only non-relativistic matter (cold dark matter and baryon), we havepm = 0. By comparing the above modified 
Friedmann equations with the ordinary ones in general relativity one can obtain the energy densities and pressure 
for f{T) gravity as [26] 



(30) 
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Figs. 8a and 8b represent the evolution of the energy density pr and pressure pr respectively against cosmic time t. 




Fig.9 



Fig. 10 



and 



Fig. 9 shows the evolution of the equation of state wt = — with cosmic time t. 
Figs. 10 represents the evolution of /(T) with cosmic time t. 



Pt = 



2k2 



-8HTf"{T) + {2T - 4:H)f{T) - f{T) + AH + QH' 



and the equation of state 



Wt = ^1 



SHTf"{T)+4Hf'{T) 
/(T) - 2T/'(r) 



(31) 



(32) 



As we are considering the emergent universe scenario, we use the scale factor a as in equation (6). Subse- 
quently, H gets the form in (7). Following reference [49] we take /(T) as 



/(T) ^jT + XT'"' 



(33) 



where, m > and T = —6H^. Consequently, 



H = — 



;3(6i/2(l + 7) - (-i72)m(6™ + 2i+"3"(-l + m))mA) 



(34) 
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The above equation is non-linear in H and hence it is difficult to get an analytical solution of a. Thus, like 

the f{R) gravity we take the scale factor in the form (6) to consider the emergent universe scenario. Using the 
above forms of a and f{T) we reconstruct pr, Pt and the equation of state parameter wt- Forms of pT and pt 
are obtained as follows 



+e3*"7 + 3e2*"/37 + 24e2*"n2a3/37) (35) 
+(3e*« + /3)/327 - 22+'"3™ma/3(e*" + PfXif}] 



_ K^r 2e'°na''(r»e'°(3+27)+2;3) ^ 1 
PT — 2^ (e'-+^i)2 + (et<»+/3)5 

+2i+"3™e2*"(-l + m)j3 + 6'"e*"(-l + Am)j3^ + 2i+'"3™m/33)A))] 

where, 



(e*«+/3)2 



(37) 



Figures 8a and 8b show the variation of pr and pr with cosmic time t. Which show that in the emergent 
scenario under /(T) gravity the energy density increases with time, while the negative pressure falls with 
increase in the cosmic time. The equation of state parameter wt < — 1 as plotted against cosmic time in figure 
9. This indicates the phantom scenario. In figure 10 we find the evolution of f{T) with cosmic time. 



VI. CONCLUDING REMARKS 

In this work, we have considered emergent universe in generalized gravity theories like Chameleon, f{R) and 
/(T) gravities. While considering the Chameleon field, we have considered flat, open and closed universes. 
We have reconstructed the potential of Chameleon field under emergent universe scenario and observed its 
increasing nature with evolution of the universe. Also, we observed the increasing behavior of the associated 
scalar field in all of the flat, open and closed universes. We have derived the conditions for emergent scenario 
in the situations of flat, open and closed universes. Next, we have considered the emergent universe scenario 
under f{R) gravity. Taking a particular form of f{R) we have considered the energy density, pressure and 
subsequently the equation of state parameter. It has been revealed that the equation of state parameter 
exhibits quintessence like behaviour (above —1). However, considering the behaviour of the equation of state 
parameter under /(T) gravity in the emergent universe scenario, we have seen it to lie below —1. This indicates 
phantom like behaviour. 
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